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ABSTRACT 

The dynamical behavior of two types of non-equilibrium systems is discussed: (a) two- 
dimensional cellular structures, and (b) living polymers. Simple models governing their 
evolution are introduced and steady state distributions (cell side in the case of cellular 
structures and length in the case of living polymers) are calculated. In both cases the 
models possess a one parameter family of steady state distributions. Selection mechanism 
by which a particular distribution is dynamically selected is discussed. 
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Pattern selection is well known to take place in systems far from thermal equilibrium 
such as growing crystals, reaction diffusion problems and in many others. In these sys- 
tems a particular pattern, or structure, is dynamically selected out of a family of possible 
structures. The canonical problem for which selection has been studied and explicitly 
demonstrated is that of a front propagation, where the velocity with which a stable state 
propagates into an unstable one is selected [1,2]. This mechanism was then applied to 
study the pattern selection occuring in systems quenched beyond their limit of stability 
[3-8]. 

In the present paper we briefly review two interesting classes of systems which have 
recently been suggested to display dynamics governed by selection. The first class includes 
systems with cellular structures [9,10] and the second is living polymers [11]. Here the 
selected quantity is not a spatial structure but rather the steady state distributions of 
certain quantities like that of the cell side in the case of cellular structures and the polymer 
length distribution in the case of living polymers. This review summarizes results obtained 
inRefs. [9-11]. 

A. Two Dimensional Cellular Structures 

Cellular structures are rather common in nature [12]. Examples include polycrystal 
[13 — 15], foams, soap froths [16,17], magnetic bubbles [18 — 21], and many others. In 
many cases these structures are metastable and they keep changing in time, evolving 
towards equilibrium. For example in the case of polycrystals, the cells are composed of 
single crystallites which are randomly oriented. Neighboring cells are separated by grain 
boundaries which cost some energy. In an attempt to reduce this energy the cells coarsen, 
thus reducing their area and the energy of the grain boundaries. Similarly, foams are 
made up of liquid film membranes separating gas filled cells. The energy of the system 
is given by the surface energy of the membranes. Here too the energy is reduced via 
the coarsening process. The main difference between the two systems is that while the 
surface free energy of liquid membranes is isotropic, the grain boundary energy strongly 
depends on the relative orientation of the neighboring crystallites and on the direction 
of the boundary. This certainly affects the details of the dynamics, but not the general 
tendency of coarsening. The evolution of these systems leads to a scale invariant state, with 
steady state distributions of some of their properties. Most studies so far have concentrated 
on two dimensional (2D) systems which are much simpler to investigate both theoretically 
and experimentally than in 3-D [22 — 29] . Here we discuss the evolution of 2D soap froth. 

Consider a two dimensional array of soap bubbles made up of liquid film membranes 
separating gas filled cells. The edges of the cells meet at vertices with coordination number 
q = 3. The reason is that any vertex with higher coordination number tends to split into 
several vertices, generically reducing the energy, or the total edge length, of the configura- 
tion. If the gas cannot diffuse through the membranes, the system reaches a mechanical 
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equilibrium in which the angle between any pair of edges meeting at a vertex is 120°. How- 
ever, due to the permeability of the membranes, the gas diffuse from one cell to another 
and the system slowly evolves in a way such that its total edge length is reduced. In 1952 
von Neumann considered the evolution of 2D soap froth under the assumption that the 
diffusion process is sufficiently slow so that the system is at a mechanical equilibrium at 
any given time. He showed that in this limit, the area of an / sided cell, Ai, satisfies the 
following equation [22] 

dAl K(l-6) (1) 



dt 

where K is a constant which is proportional to the surface tension and to the permeability 
of the membranes. This is a very important and surprising result, since it shows that the 
evolution of the area of a cell is independent of the details of its neighborhood (such as the 
area, the shape and the number of edges of the adjacent cells), but only on the number 
of edges of the cell itself. Eq. (1) suggests that / < 6 sided cells shrink and eventually 
disappear while cells with / > 6 sides grow in time. The area of hexagonal cells remains 
unchanged. As a result of the disappearance of an / < 6 cell, the number of sides of its 
neighboring cells change (see Fig. la). When the disappearing cell is a triangle, each of its 
adjacent cells loose a side. For a disappearing square, two of its neighbors loose a side 
and two remain with the same number of sides. In the case of a pentagon, two neighbors 
loose a side, one gain a side and two remain unchanged. These events are known as T2 
processes. As a result of these processes, / < 6 sided cells do not disappear from the system 
altogether. This is consistent with Euler's law which states that V — E + F = 2 where V is 
the number of vertices, E is the number of edges and F is the number of faces of a closed 
two-dimensional graph in a plane. Using the fact that in our case three edges meet at each 
vertex, it is easy to show that the the average number of sides per cell is 6. Therefore, any 
configuration in which / > 6 sided cells are present, must include / < 6 sided cell too. In 
addition to the T2 processes, one also finds Tl processes corresponding to edge switching 
events (see Fig. lb). In soap froth these are rather rare events and therefore we will not 
consider them here. Also, two sided cells are metastable and they are not generated during 
the evolution. We therefore assume that the system is composed of / > 3 sided cells only. 



(a) 



Tl 



(b) 



Figure 1: (a) T2 processes, describing the disappearance of 3, 4 and 5 sided cells, 
respectively. (b) side switching Tl process. 

To study the evolution of soap froth we consider Ni(t), the number of /- cells (i.e. I- 
sided cells) at time t. Let iV = Ni be the total number of cells in the system, and 

define 

x^Nt/N (2) 

Clearly, the concentration x\ of I- cells satisfy 

00 

= 1 (3a) 

1=3 

Moreover, due to Euler's law one also has 

00 

^Zx z =6 (36) 

1=3 
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To proceed, we introduce three parameters 103, w 4 and w$ which yield the vanishing rates of 
triangles, squares and pentagons, respectively. These rates together with the T2 processes 
determine the evolution equations of Ni. For example the resulting equation for N 3 takes 
the form 

^ = -w 3 N 3 + [3P fl (4, S)w 3 N 3 + 2P fl (4, A)w 4 N 4 + 2P fl (4, 5)w 5 N 5 ] - Pa(3, 5)w 5 N 5 (4) 

where Pn{l,m) is the probability that an edge of a disappearing m < 5-cell is shared 
by an adjacent /-cell which undergoes a side removal process. Similarly, Pa(1,5) is the 
probability that an edge of a disappearing pentagon is shared by an adjacent /- cell which 
undergoes a side addition process. In this equation the first term represents the rate at 
which triangles disappear due to the von-Neumann evolution. The next three terms yield 
the rate at which squares turn into triangles due to the disappearance of a neighboring 
3, 4 and 5 cells, respectively. The last term gives the rate at which triangles gain a side 
and become squares due to the disappearance of a neighboring pentagon. 

We now introduce the mean field approximation for Pr(1, m) and Pa(1, 5). In consid- 
ering Pa(1, 5) we neglect all correlations between neighboring cells and take 

P A (Z,5) = ^ = ^ (5) 
v ; 6N 6 y J 

In Pn(l,m) some correlations have to be retained. In particular, in order not to generate 
2-cells a side removal process should not take place in a triangle. We therefore take 

Pi?(Z ' m) = {/V(6-3x3) i />3 (6) 

where the denominator (6 — 3x3) is taken to ensure that the normalization condition 
YlhL 3 PrQ"> m) = 1 is satisfied. One can easily write down similar equations for Ni, I > 3. 
Summing these equations one obtains an equation for the decreasing rate of the total 
number of cells 

dN 

— = -(w 3 N 3 + w 4 N 4 + w 5 N 5 ) (7) 

Combining this equation with the equations for Ni one obtains the following equations for 

xi, I > 3 

dx 1 

— = -w t xi + K xi + Kx[(l + l)xi +1 - ailxi] - K 2 [lx t - a t (l - (8) 

where 

K = w 3 x 3 + w 4 x 4 + w 5 x 5 

K\ = (3w 3 x 3 + 2w 4 x 4 + 2w 5 x 5 )/(6 - 3x 3 ) (9) 
K 2 = w 5 x 5 /6 
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and 

f , 1=3 
ai= \l , l>3 

Also -u;; = for Z > 5. 

We are interested in the steady state solutions of Eqs. (8,9) for which dxi/dt = 0. It 
is easy to demonstrate that these equations have a one parameter family of steady state 
solutions. To see that this indeed is the case we notice that for any given < xs < 1 the 
first two equations (corresponding to / = 3, 4) could be solved to yield X4 and x$. One 
can then solve iteratively for the rest of the distribution, finding xi from the (/ — l)-th 
equation. The family of distributions obtained in this way may thus be parametrized by 
x 3 . 

We now consider the steady state distributions in more detail. Let 

VI = Ixi (10) 
The steady state equations for / > 5 take the form 

K lVl+1 - (K, + K 2 - ^) yi + K 2Vl _ x = (11) 

These equations have solutions of the form 

yt = \ l l K (l + 0(j)) (12) 

Inserting this form into Eq. (11) the parameters A and k may be obtained. Two solutions 
are found. One with 

K 2 K 
Al = K? K1 = K^W 2 (13a) 

and the other with 

The general solution, therefore, takes the form 

Xl = A(J^J l^-i + Bl"*- 1 (14) 

The full distribution (including the / < 5 concentrations) is found by first chosing a par- 
ticular value for X3 and solving the / = 3,4 equations to obtain x^ and £5. One then uses 
the form (14) for the / > 5 cell concentrations and determine the parameters A and B by 
demanding that the two sum rules (3a and 36) are satisfied. The tail of the distribution is 
thus given as a sum of two terms: an A term which decays exponentially with / and a B 
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term which decays algebraically. By varying x 3 the whole family of possible steady state 
distributions is obtained. 

In calculating these distributions it is found that B is positive for sufficiently small X3, 
and it changes sign at some point X3 = x$. The distributions with B < are unphysical 
since xi become negative for large /. Moreover, linear stability analysis shows that those 
solutions with B < are unstable while those with B > are stable. As in other problems 
where selection takes place it is expected that the physical solution which exhibits the 
fastest decay of x\ with / is the one which is selected. Here it is the B = solution, which 
is the only one exhibiting exponential decay. All other physically relevant distributions 
(with B > 0) display algebraic decay for large /. Numerical studies of the dynamical 
equations indicate that indeed this is the case and that the system evolves towards the 
B = steady state distribution. The resulting distribution agrees well with experimental 
results obtained in soap froth [9, 10]. 

In order to check how good is the mean field approximation applied in this work, we 
consider now the case in which only Tl processes take place. Here it is assumed that the 
gas does not diffuse from one cell to another, hence T2 processes do not take place and 
cells do not disappear. Thus, the number of cells in the system is conserved throughout 
the evolution. This problem has been solved exactly [30] , yielding a distribution which for 
large / takes the form 

xi ~ X l /Vl (15) 

with A = 3/4. 

To model this process within the mean field approximation, we define the rate w as 
the probability that a Tl process takes place per time unit per edge. Each such event 
results in a removal of a side from two cells and addition of a side to two others. The 
process cannot take place if either of the two neighboring cells which lose an edge in this 
process is a three-cell. This ensures that no two-cells are generated in this dynamics. The 
evolution equations of x\ are given by 

^ = w^^[(l + l)x l+1 - l Xl ] + w (^^) 2 [(/ - - lx t ] (16) 

These equations process a single fixed point distribution 

xi ~ X l /l (17) 
with A ?a 0.85. This qualitatively agree with te exact answer. 

B. Living Polymers 

Another example of a dynamical system whose steady state distribution is uniquely 
selected out of a family of distributions is given by living polymers. Here we consider a 
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system of aggregates (say, polymers or micells), each composed of n units. Let c n be the 
number of n- mers in the system. Under equilibrium condition the system reaches a well 
defined length distribution c n . We consider the case in which the system is driven out 
of equilibrium by pulling out monomers (or small size aggregates) at some rate. Clearly, 
c n decrease with time and the system eventually vanishes. The question is how does the 
length distribution c n look like in the long time limit. To study this problem we consider 
the concentrations x n = c n /c where 

oo 

c=J^nc n (18) 

n=l 

and ask whether these quantities reach a well defined limit for large time. 

As in the case of cellular structures, one finds that the dynamical equations which 
govern the evolution of x n possess not just a single fixed point but rather a one parameter 
family of steady state distributions. Again a dynamical selection takes place by which a 
particular distribution is selected. To study this behavior we introduce a model describing 
the evolution of living polymers. For simplicity we assume that the interaction between the 
various n-mers take place via association-dissociation process in which an n-mer absorbs 
or emits a single monomer and becomes (n + 1) or (n — l)-mer, respectively. The equations 
governing this process take the form: 

d oo oo 

—A = -kxi ^2c n + k^c n - aci, 

n=l n=2 

^ = fcxi( Cl /2 - c 2 ) + k(c 3 - c 2 /2), ( 19 ) 
dc n 

—j^- = kxi{c n -\ - c n ) + k(c n+1 - c n ) for n > 3. 

where k is the dissociation rate, kx\ is the association rate of an n-mer with a monomer 
and a is the rate at which monomers are pulled out of the system, x\. The association 
rate is taken to be proportional to the concentration of the monomers in the system. For 
simplicity it is assumed that k and k are independent of n. This assumption does not 
change the qualitative results derived here. The equation for c 2 is somewhat different from 
equations for c n with n > 2. First, the rate at which two monomers combine to yield a 
dimer is kxi/2 and not kx\. In addition, while a dissociation of n-mer with n > 2 can 
take place near either of its two ends, in the case of a dimer the dissociation can take place 
only at one point. Hence the factor of a 1/2 in its dissociation rate. Summing Eqs.(19) 
one finds 

dc 

- = -ac u (20) 
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This equation yields the rate at which the number of units in the system decreases with 
time. 

Combining Eqs.(20) with (19) one obtains the following equations for x n . 

d °° 00 

= —kX\ x n + k X n - ax l + ax l, 

n=l n=2 

= kx 1 (xi/2 - x 2 ) + k(x 3 - x 2 /2) + ax 1 x 2 , 

= kxi(x n -i - x n ) + k(x n +i - x n ) + ax\x n for n > 3. 

This is a set of non-linear equations with non-local interactions of the type found in the 
equations describing the evolution of cellular structures. Here the "head" of the distri- 
bution, xi, directly interacts with its "tail", x n , for arbitrarily large n. The fixed point 
distribution of these equations yield the long time behavior of c n . Note that the normal- 
ization 

00 

J2nx n = l. (22) 

n=l 

is preserved by Eqs. (21). 

Eqs. (21) possess a one parameter family of fixed point distributions. To see that this 
is the case one first makes an arbitrary choice for x\ and x 2 . The second equation of (21) 
determines X3. One can then successively solve the remaining equations, obtaining x n for 
n > 3. The two free parameters x\ and x 2 may now be varied so as to satisfy the equation 
for x\ (or equivalently, the normalization condition (22)), leaving one free parameter, say, 

X X . 

The fixed point distributions may easily be obtained by noting that for a given xi, 
Eqs. (21) are linear in x n , n > 2. Setting dx n /dt = one finds that the fixed point 
distributions are of the form 

x n = A\ n _~ 2 + B\ n + ~ 2 , n > 2, (23) 
where A± are the roots of the characteristic equation 

k\ 2 — (k + kxi - axi)X + kx\ = 0. (24) 
Here A and B are parameters which are determined by the first two equations in (22) 

A B _ (k - a)x\ + ax\ 



1 — A_ 1 — A+ k — kx\ 

(25) 

A\k\- - (k/2 + kx x - ax-C)] + B[k\ + - (k/2 + kx x - ax x )\ = -kx\/2 . 
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Eqs. (25) yield the two parameters A and £?, thus determining the length distribution. 
Clearly, all x n 's obtained in this way must be non- negative for the distribution to be 
physically meaningful. For this to be the case one has to require that (a) the roots X± 
are real with < A_ < A + < 1, (b) X2 = A + B > and (c) B > 0. These requirements 
ensure that x n decays with n purely exponentially without oscillations. When the decay is 
oscillatory some of the x n 's become negative and thus unphysical. The B > requirement 
is needed since B corresponds to the slower decay rate of x n . A negative B yields negative 
x^s for large n. 

We now consider the question of selection of a unique steady state distribution out 
of this family. By examining Eq. (24) one can verify that for any set of parameters k, k 
and a there exists an xm such that as long as < x\ < %, condition (a) is satisfied. 
It is also easy to verify by solving Eq. (25) that condition (b) is satisfied as long as 
< x\ < xm- As for condition (c), one can show that the parameter B is positive for 
small x\ and changes sign at x\ = xs, where xs depends on the parameters k,k and a. 
For xs > xm all fixed points with < xi < xm are physically relevant. On the other hand 
for xs < xm only those with < x\ < xs correspond to physical distributions. These 
two cases yield two different selection mechanisms. Following the work of Aronson and 
Weinberger [1], and Dee et al [3] one may conjecture that in the first case the marginal 
fixed point corresponding to x\ = xm is selected, in the sense that as long as the initial 
length distribution decays sufficiently fast with n, the system evolves towards the xi = xm 
fixed-point. On the other hand, in the second case (the non-linear marginal stability, or 
case II) the selected fixed-point is the one corresponding to xi = xs- In both cases the 
selected distribution corresponds to the physically accessible fixed-point with the fastest 
decay rate of x n with n. 

This model has been studied numerically, and it was shown that indeed the long time 
behavior of the system is governed by either the x\ = xm or the x\ = xs fixed points 
depending the parameters k, k and a. One of the parameters, say k, may be taken as 1. It 
determines the time scale in the problem and does not affect the steady state distribution. 
The resulting (/c, a) phase diagram is given in Fig. 2. 
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Figure 2: The (k,a) phase diagram of the model (21) for k = 1. Two regions are 

found. One (M) in which the marginal fixed point x\ = xm is selected, 
and the other, (S), in which x\ = xs is selected. 
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